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Abstract UVOIR bolometric light curves provide a valuable insight into the nature of type la supernovae. We present an 
analysis of sixteen well-observed SNe la. Constraints are placed on several global parameters concerning the progenitor system, 
explosion mechanism and subsequent radiation transport. By fitting a radioactive decay energy deposition function to the quasi- 
exponential phase (50 to 100 days after maximum light), it is found that the ejected mass varies by at least a factor of two. This 
result suggests that a sub Chandrasekhar mass white dwarf could be the progenitor system of some SNe la. We find that the 
range in the amount of synthesized 56 Ni indicates a significant variation in the burning mechanism. In order to explain a factor 
of ten range in the observed bolometric luminosity, more detailed modeling of the explosion mechanism is required. 
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Type la supernovae (hereafter SNe la) have become an excep- 
tional tool in modern cosmology. Due to their large luminos- 
ity, they are used to place constraints on cosmological param- 
eters, and as of yet, provid e the only direct evidence for the 
existence of dark energ y jRiess et alJll998at IPerlmutter et alJ 
: ll999llLeibundgutl200ll) . Despite the insight SNe la have given 
1 us about the universe, several key issues related to the nature 
of their progenitor system(s) and the physics of the explo- 
sion mechanism(s) have remained unsolv ed (for reviews see 
iHillebrandt & NiemeveJ200filLiviol2000l) . 

Today it is commonly believed that SNe la emerge from 
the thermonuclear incineration of a carbon oxygen (C-O) 
whit e dwarf e xploding near or at the Chandrase khar mass 
jNomoto . Thielemann. & Yokoil 1 19841 IWooslev & WeaveJ 
1986; Hillebr andt & NiemeveJ EoOO). The energy released 
from burning to nuclear statistical equilibrium (NSE) com- 
pletely destroys the white dwarf. The optical/IR light curves 
are powered by the Comptonization of 7 rays produc ed from 
the radioactive decay chain 56 N i — > 56 Co — ► 56 Fe JPankevI 
1 19621 IColgate & McKeel \l969). Within this paradigm the 
C-O white dwarf accretes matter from an evolved massive 
star either via Roche lobe overflow or through stellar winds. 
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With the detection of Ha in the well-observed SN 20 02ic 
jHamuv et all2003r.lNomoto et all2004t.lKotak et all2004l) the 
single degenerate model has become favored over the double 
degenerate model. 

It is clear that SNe la are not a homogeneous group of 
stellar explosions rather, they display a range in luminosity of 
a factor of ten or more ( Suntzeffj fl996t IContardo etaTl feoOO: 
Suntzeff 2003; Stritzinger 2005). In fact not one self-consistent 
explosion model has yet been presented that can successfully 
account for this observed range in luminosity. This lack of un- 
derstanding of the progenitor systems is unsettling, and must 
be addressed, if we are to have confidence in the cosmological 
results provided by SNe la. 

Fortunately, during the past de cade, a number of observing 
campaigns (see Leibundgut 2000, for a list) have obtained ex- 
cellent data sets for a large number of SNe la. With these data 
sets we are currently in a position to conduct a systematic in- 
vestigation of their photometric and spectroscopic properties. 
The motivation of this work is to place further constraints on 
both the progenitor systems and explosion mechanisms with 
these observations. In particular, we construct UVOIR bolo- 
metric light curves from broad-band photometry of a number 
of SNe la. With the UVOIR light curve we determine several 
global parameters, e.g. total ejected mass, 56 Ni mass and the 7- 
ray escape fraction. Investigating the range of such parameters 
can provide insights into the nature of SNe la, and help differ- 
entiate between the various paths of stellar evolution that the 
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progenitors follow, as well as the manner in which thermonu- 
clear combustion occurs. 

The structure of this paper is as follows. Sect. 2 discusses 
the observational data used and the UVOIR light curves. Sect. 3 
describes our analysis of the UVOIR light curves. Sect. 4 con- 
tains the results of this analysis and we conclude in Sect. 5. 

2. Observational data and UVOIR light curves 

(U)BVRI-band observations for sixteen SN la have been pro- 
cured from a variety of sources. Table 1 lists all events consid- 
ered in this study along with references to the sources of the 
data, as well as a number of important parameters used and 
determined in this work. All events listed in Table 1 have ex- 
cellent optical photometry that extends fro m pre-max imum out 
to ~ 100 days past maximum light. (See IStritzingerlEo05l for 
more details concerning the method and sample.) 

To determine the amount of 56 Ni produced in the ex- 
plosion from the observed UVOIR flux, two parameters 
are required. They include an estimate for the total extinc- 
tion and the distance to the host galaxy. Values listed for 
G alactic r eddening are those given by the COBE dust maps 
of Schleg eletail 11 998). Host galaxy reddenings were selected 
from Phi llips et alJ ( 119991) for those SNe la that coincide with 
our sample. For more recent events we adopted host galaxy ex- 
tinctions from the literature, with preference to those calculated 
by the Phillips method. The distances chosen depended mainly 
on what was available in the literature. For SNe la without a 
direct distance measurement, in this case, a Cepheid distance 
measurement or a surface brightness fluctuation (SBF) distance 
we used a Heliocentric velocity obtained from NED and con- 
verted this to the CMB reference frame. 1 For each CMB dis- 
tance we have assumed an uncertainty of 300 km s _1 . 

With the advent of the Hubble Space Telescope (HST) there 
has been a substantial effort from two groups (namely the HST 
Key Project (HKP) and the Saha, Tammann and Sandage (STS) 
group) to obtain accurate Cepheid distance measurements to 
galaxies that have hosted SNe la. Typically these two indepen- 
dent groups determine different distances for any one galaxy, 
even though they use the same data and similar data reduc- 
tion software. These differences are a reflection of assumptions 
made in their analysis. The main factor that contributes to these 
discrepancies is the exact P-L relation used, and several other 
subtle nuances exacerbate the problem. These include: (1) the 
criteria adopted to select the Cepheids used to determine the 
distance, (2) if and how metallicity corrections are applie d, and 
(3) anomalies related to the camera(s) on H ST (see lParodi et alJ 
2000: lGibson et all200(AlRiess et all20()3) . 

In short, the HKP obtains a short distance scale that 
leads to a value of H ~ 70 km s -1 Mpc _1 while the STS 
group determine s a long distance scale that gives H ~ 60 
km s _1 Mpc -1 . Riess et al. (2005) have reviewed this issue, 
in order to reconcile the distances obtained by the two groups 
to galaxies that have hosted a SN la. 

1 Through out this work we have adopted 
H = 72kms" 1 Mpc" 1 . 



Four of the SNe la used in our study have a direct Cepheid 
distance to their host galaxy and one has a Cepheid distance 
to its galaxy cluster. Whether the distance scale is long or short 
will depend on which data (i.e. from which of these two groups) 
we use. This in turn will lead to either an under or over estimate 
of the distance. This uncertainty in the distance will affect the 
a6 Ni mass we determine and hence the ejected mass. 

To illustrate the effect of this on our results, we consider the 
galaxy NGC 3982. This galaxy has three independent Cepheid 
distance measurements, one from each of t he two teams just 
mentioned and one from iRiess et all d2005l) who have used a 
new calibration of the P-L relationship and an elegant metal- 
licity correction. 

The STS group has publis hed a d istance modulus of 
p = 31.72 ± 0.14 ISahaet all 12001 bb which is based on 
the Ce pheid P-L relation published by iMadore & Freedmanl 
11991b . Using a P-L relation based on ^650 Cepheids lo- 
cated in the Large Magellanic Cloud that were observed 
by the Optical Gr avitational Lensing Experiment (OGLE) 
llJdalski et alJ 1 1999b the HKP t eam concluded that u = 
31.56 ± 0.08 IStetson et alJl200lb . More recently iRiess et alJ 
(2005) used the Advanced Camera for Surveys on the HST to 
obtain a Cepheid distance. They use d th e P-L relation pre - 
sented by Tammann & Reindl (2002) and iThim et alJ (|2003), 
based on only those Cepheids in the OGLE data set that have 
periods longer than 10 days. With this new P-L re lation, and 
after applying a metallicity correction, IRiess et all found p = 
31.66 ±0.09. 

With the HKP distance modulus we obtain a 56 Ni mass that 
is 9% le ss com pared to the STS distance modulus, and 4% less 
with the lRiess et al J distance modulus. In turn the ejected mass 
derived with the HKP distance modulus is 22% greater com- 
pared to_the_STS distance modulus, and 13% greater compared 
to the lRiess et al J distance modulus. 

We have used a SBF distance for three SNe la in this sam- 
ple. The zero-point for this method is based on an empirical 
relation derived from Cepheids, which can lead to a systematic 
under- or over- estimation of the distance, depending on whose 
Cepheids are used for the calibration. The three SBF distances 
w e have used are calibrated with the HKP Cepheid distances 
of lFerrarese et alJl2000l) . As pointed out bv lTonrv et all 1200 lb 
the zero-point for the SBF distance scale is still being im- 
proved. With the application of the Cepheids determined with 
the complete set of OGLE data the SBF distances used here 
would be ~0. 1 mag fainter. With this correction our derived 
56 Ni masses would increase by ^8% and the ejected masses 
would decrease by ~16%. 

Finally, in Table 1 we list two parameters calculated in this 
study: (1) an estimate of the amount of 56 Ni produced from 
burning to NSE (see Sect. 13. U . and (2) an estimate of the time 
when the ejecta makes the transition from being optically thick 
to optically thin (see Sect. 13. 2> . 

To construct the UV OIR light curves we u s ed the metho d 
already employed by Vacca & Leibundgut ( 1996, 1997), 
IContardo et alJ 12000b and Stritzing er & Leibundgut! 12005b . 
We refer the reader to these articles for a more detailed descrip- 
tion of this empirical fitting method, which we briefly summa- 
rize here. 
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Each filtered light curve is fitted with a ten-parameter func- 
tion. This function consists of a Gaussian for the peak phase, 
a linear decline for the late-time decay (i.e. 56 Co — > 56 Fe), an 
exponentially rising function to fit the initial rise to maximum, 
and a second Gaussian for the inflection or secondary maxi- 
mum that is observed in the VRI-bmd light curves. 

SNe la without a £/-band light curve, a correction was 
added, as described by IContardo et alJ J2000I). IC ontardo etalj 
used a correction based on SN 1994D jRichmond et al.H l 995: 
iMeikle et al.lll996llSmith et allEo OO). however, this event had 
an unusual blue color at maximum. For this reason, correc- 
tions based on SN 1994D tend to overestim ate the fraction o f 
flux associated with the t/-band photometry ( Stritzinger 2005). 
Inste ad we em ployed a correction derived from SN 1992 A 
JSuntzefJ Il996l) . Although there are many well observed 
events, SN 1992 A is one of the only normal SN la that does 
not suffer extinction due to host galaxy reddening. 

Owing to a lack of data we did not include the fraction of 
flux associated with wavelengths above 10,000 A. However, 
at maxi mum light the JH A'-bands contribute no more than 
~ 5% JSuntzefl 19961) to the total bolometric flux. Around 60 
days after maximum light, when the bolometric light curve fol- 
lows nearly a linear decline, the infrared contributio n rises to 
no m ore than ~ 10% of the total bolometric flux JCon tardo 
1200 ll) . 

To produce the UVOIR light curve, each fitted light curve 
is converted to flux (ergs s _1 ). Next a reddening correction 
is applied, and then each filtered light curve is summed to 
obtain the total flux. Note that we do n ot normalize the flux 
to any decline rate relatio n (e.g. Ami ^ ( Phillips et al.l| 1 999l). 
MLCS jRiessetal.|[l996l) or stretch dPerlmutter et all 1 19971 
iNobili et all2005l) . 

3. Determining global parameters of SNe la 

In this section we describe the manner in which UVOIR light 
curves are used to determine the parameters of interest here. 

3. 1. 56 Ni mass 

At maximum light the peak luminosity of a SN la is related, to 
first order, to the amount of 56 Ni produced during the explo- 
sion. The amount of 56 Ni synthesized from burning to NSE is 
itself thought to be largely dependent on the explosion mech- 
anism. With the 56 Ni mass we are directly probing the most 
sensitive part of the explosion and can use observations to place 
const raints on the e xplosion mechanism. 

ISuntzefJ Jl996l) showed that at maximum light ^80% or 
more of the total flux from a SN la is emitted in the op- 
tical. Therefore with UVOIR light curves, constructed from 
UBVRI broad-band photometry, one easily obtains a measure 
of the total flux and, through application of Arnett's Rule, the 
56 Ni mass. Arnett's Rule states that at maximum light, the lu- 
minosity of a SN la is equal to the instantaneous energy de- 
position ra te from the radioactive de cays within the expand- 

I ll i I II I 

ing ejecta llArnettl ll982: Arnett et al. 1985). To determine the 

56 Ni mass we use the simple relation that gives for 1 Mq of 



56 Ni, a total luminosity at maximum light of 

L max = (2.0 ± 0.3) x 10 43 -P- ergs" 1 . (1) 

Mq 

The error in Eq. Q corresponds to a 3-day uncertainty in 
the a dopted b olometric rise tim e of 19 days (see Sect. 4.1 of 
IStritzinger & Leibundeutl2005l for more details). As approxi- 
mately 10% of the total flux at maximum light is emitted out- 
side of the optical, each 56 Ni mass derived from Eq. (1) has 
been increased by a factor of 1.1. The dominant errors in the 
deduced 56 Ni mass are associated with the adopted distance to 
the host galaxy and the total extinction dContardo et al.lZo'ooT) . 

3.2. Total ejected mass 

To place constraints on the ejected mass we perform a least- 
squares fit of a radioactive beta-decay energy (RDE) depo- 
sition function to the post maximum phase UVOIR light 
curve. Prior works that discuss this method inc lude the pio- 
neering investigations of IColgate et al. Nl980alb | ). followed by 
the more sophistic ated treatment presented by Jeffery ( 1999); 
(see also ICappellaro et al.Nl997| iMilne. The. & Leising| ll999 
and lMilne. The. &L eising 2001 for similar methods and tech- 
niques.) However, as of yet, no attempt has been made to apply 
such a method to UVOIR light curves derived from real obser- 
vations. 

An expression for the energy deposition of N^io atoms of 
56 Ni in the optically thin limit (i.e. when r 1) is represented 

by 



Ed C p = E m + E Co e + + [1 - exp(-r)]A Co 7 
= AniNnio exp(— ANit)Qisri 7 

+ AcoNmot — ~~~t — [exp(-Acot) - exp(-A Ni t)] 

x {<3coo+ +Qco 7 [l -exp(-r)]} . (2) 

Ani and Ac are the respective e-folding decay times of 8.8 
and 111.3 days for 56 Ni and 56 Co. Q Ni 7 (1.75 MeV) is the 
energy released per 56 Ni — > 56 Co decay. Q Co c + (0.12 MeV) 
and Qco j (3.61 MeV) are the positron and 7-ray energies re- 
leased per 56 Co — ► 56 Fe decay (for a detailed discussion of the 
radio active properties of this decay chain see, e.g. Nadvozhin 
Il994h Note that throughout this work we assume that all neu- 
trinos produced from the 56 Ni^ 56 Co^ 56 Fc decay chain es- 
cape the ejecta entirely and do not contribute to the observed 
UVOIR flux. 

As Eq. (I3 is only applicable in the optically thin limit, 
when the thermalized photons can freely escape, it is safe to 
assume that at these epochs the majority of 56 Ni has decayed 
to 56 Co, and therefore the remaining amounts of 56 Ni provides 
a negligible contribution to the energy deposition. At these 
epochs the UVOIR light curve appears to be nearly first order 
exponential, howe ver it is more accurately described as 'quasi- 
exponential' (seeljeffery 1999, for a detailed discussion). With 



4 



M. Stritzinger et al.: Bolometric properties of SNe la 



the presence of only one radioactive species, the mean optical 
depth r has a simple t~ 2 dependence: 



(3) 



If we replace r in Eq. (|2ji with Eq. Q, and then perform a 
least-squares fit of Eq. to the UVOIR light curve (between 
50 and 100 days past maximum light when Eq. Q is valid), we 
can determine the 'fiducial time' t . It is at this time that the 
ejecta becomes optically thin. 

Following the discussion o fl.Teffervi ri99qi. one finds that t 
can be expressed as 



MejKq 

8tt 



(4) 



The variable M C j is the total ejected mass, n is the 7-ray mean 
opacity, v e is the e-folding velocity of an exponential model's 
density profile, and q is a general form factor that describes the 
distribution of 56 Ni in the ejecta. 

During the optically thin phase for an all-metal ejecta (/i e = 
2), k is expected to be in the rang e 0.025 to 0.033 cm 2 g _1 (see 
LSwartz et alJfl995t |Jefer3[l999, and references therein for a 
detailed discussion). We adopted the value of 0.025 cm 2 g~ x 
as our fiducial 7-r ay mean opacity. 

I-Teffervl ill 9991) lists the model e-folding velocity of sev- 
eral successful 1-D explosion models consisting of 1.4 Mq 
Chandrasekhar-size white dw arfs. These e-folding velocities 
are ~ 2700 km s _1 for W7 j^mojg^Tjiielemann. & Yokoi 
1 1984 . 2750 km s" 1 fo r DD4 JWooslev & Weaverll994al). and 
3000 k m s" 1 for M36 faoflichl 19951) In additionljefferv et all 
i ll 9921) found that the DD2 model of lWooslevI dl99lh has a e- 
folding velocity o f ~3160 km s" 1 . 

More recently iRooke & Hillebrandtl |2004) published two 
full-star 3-D explosion models of a 1.4 Mq white dwarf; with 
different ignition conditions: a centrally ignited configuration 
(c3_47r) and a foamy multi-bubble flame structure (/l). Using 
Eq. (A 10) of Jeff er^ (1 19991) and parameters given in Table 1 of 
Ropke & Hillebrand7 T2004). we have calculated the e-folding 
velocities for these two models. The e-folding velocities corre- 
spond to ~1611 km s _1 for the c3_47r simulation and ^1842 
km s _1 for the fl simulation. These values are substantially 
smaller than the previously cited 1 -D models and reflect the dif- 
ference between the density profiles generated by 1-D and 3-D 
simulations. Note that all these models are based on the explo- 
sion of a 1.4 Mq Chandrasekhar-size white dwarf. In the cal- 
culations presented below we arbitrarily adopted 3000 km s _1 
as our 'average' fiducial e-folding velocity. 

The parameter q is equal to one for high concentrations of 
56 Ni at the center of the ejecta, small for low concentrations 
within the center, and one-third for the case when the 56 Ni is 
evenly distributed throughout the ejecta (see Jefferv 1999, for a 
detailed discussion). There is mounting evidence that an appre- 
ciable amount of 56 Ni is moderately mixed within the ejecta. 
However, it is likely that the amount of mixing may vary sig- 
nificantly from supernova to supernova. 

An analysis of early-time spectra of SN 199 IT 
jRuiz-Lapuente et alii 19921 Bylazzali et alJll995|) indicates the 



existen ce of an outer shell of 56 Ni. In contrast, iGeorgii et alJ 
(2002) presented observations of SN 1998bu obtained with 
COMPTEL. They concluded that their non-detection of 7 rays 
from the 56 Co— > 56 Fe decay chain indicates that there is no ap- 
preciable mixing of radioactive nuclides within t he ejecta in the 
context of current models. More recently Istehle et all d2005h 
have presented "abundance tomography" of SN 2002bo. With 
their unique technique they determined that the vast majority 
of 56 Ni was distributed betw een 3000 to 11,000 km s _1 for 
this particular event. LTeffervl (0j99) showed that for W7, the 
parameter q was equal to approximately one-third. As W7 
has been able to quite successfully fit observed spectra for 
normal to brig ht SNe la dHarknesslll99ll iMazzali etaill 19951 
Mazzali 2001), we have adopted a q value of one-third in the 
calculations presented below. 

With values of t , derived from the least squares fit of 
Eq. dU to the UVOIR light curve during the quasi exponen- 
tial phase, along with the adopted fiducial values for all the 
parameters in Eq. @, we can proceed to place constraints on 
the ejected mass for all of the SN la in our sample. 

3.3. j -ray escape fraction 

By comparing the UVOIR light curve to the enrgy input from 
the radioactive decays -for both cases of complete trapping of 
7 rays and complete escape of 7 rays- we can obtain a quan- 
titative description of the 7-ray escape fraction. An expression 
for the UVOIR light curve based on this prescription can be 
written as 

LC(t) ohs = (1 - 7 {t))LC(t) T>>1 + 7 (t)LC(t) T<<1 . (5) 

In this expression LC(t) b s is the UVOIR light curve, LC(t) T> i 
represents the energy input from the radioactive decays assum- 
ing complete trapping of 7 rays, LC(t) T <i represents the case 
of complete escape of 7 rays, and 7 (t) is the 7-ray escape frac- 
tion. Solving Eq. (5) for 7(f) we obtain the 7-ray escape frac- 
tion 



l(t) = 



LC(t) T >i — LC(t) b s 
LC(t) T>>1 - LC{t) 



(6) 



4. Results 



In Fig. \l\we present the least squares fits of Eq. to sev- 
eral UVOIR light curves. The four events shown in Fig. [l] are 
representative of the complete population of SN la, ranging 
from the bright SN 1991T to the subluminous SN 1991bg. Also 
plotted are the energy deposition curves corresponding to the 
56 Ni — > 56 Co — > 56 Fe decay chain for the cases of complete 
7-ray trapping (dash-dotted line) and complete 7-ray escape 
(dashed line). Table 1 lists the 56 Ni mass calculated for each 
event through Eq. Q, as well as the determined values of t . 
For this sample of SNe la, the 56 Ni mass varies by a factor of 
~10, while t varies by a factor of 1.6. 

In order to give the reader a more intuitive feeling of how 
the RDE deposition curve depends on the value of t Q , we 
present Fig. [2] This figure contains the UVOIR light curve of 
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SN 2003du, along with the energy deposition curves for dif- 
ferent values of t that vary, from top to bottom: oo, 45, 40, 
35, 32.16, 25, 20, 15, and days. As expected for a fixed 
56 Ni mass, when t is increased, the energy RDE deposition 
function evolves more slowly with respect to time. Physically 
this effect is associated with an increase in the diffusion time 
of the photons trapped within the ejecta. 

In Fig. |3 we plot t versus Ami5(UVOIR). 2 Values for 
Am 15 (UVOIR) have been determined from the UVOIR light 
curves. From this figure it is clear that there exists a corre- 
lation between these two parameters. This correlation is in 
accord with our expectations, as it is well established that 
more luminous SNe la have smaller decline rates, and thus the 
epoch in which their ejecta trans form to the ne bular phase oc- 
curs at a later time (see lPinto & Eastmanll200l|. and references 
within for a detailed discussion of the physics that describes the 
luminosity-width relation). 

Armed with our values of t , we can now proceed to place 
constraints on the total ejected mass. Fig.|4]is a plot of our cal- 
culated ejected mass versus the 56 Ni mass. For the calculation 
of the ejected mass we have used q = 1/3, v e = 3000 km s _1 
and k = 0.025 cm 2 g _1 . 

The error bars that accompany each 56 Ni mass account for 
uncertainties in host galaxy reddening and the adopted distance 
(see Table 1). For events with a CMB distance we have as- 
sumed 300 km s _1 uncertainty for (random) peculiar veloci- 
ties. 

The ejected mass error bars include: (1) the uncertainty 
listed in Table 1 for each value of t , (2) a 300 km s _1 i.e. 
10% uncertainty in v e , (3) a 10% uncertainty in k and (4) a 
30% uncertainty in the adopted value of q. These '1-er' error 
bars are not statistical but rather a sensible estimation of the 
possible range of each parameter. 

Fig. |4] displays several striking features that are worthy of 
comment. First, this figure suggests that there exists a range 
in the ejected mass of about a factor of two. Three events 
(SN 1992 A, SN 1994D, and SN 2000cx) that have moderate 
amounts of M ej (i.e. 0.4 - 0.6 Mq) are of particular inter- 
est. These events are located nearly 3-a below the most mas- 
sive events, which lie near the canonical value of 1.4 Mq. In 
order to increase the ejected mass of these three events to a 
Chandrasekhar mass, it is necessary to reduce either q (which 
is highly unlikely) or our fiducial value of k by a factor of two, 
or increase either the value of t Q by a factor of ~1.3 or v e by 
a factor of ~1.4 or more. Implementing any of these changes 
results in ejected masses for all the other 'normal' SNe la to 
be comparable to that of a neutron star mass. In other words, if 
we change any one of the parameters in Eq @ while keeping 
all others constant, there will always exist a relative difference 
in the ejected mass of ~2 between these three events shown in 
Fig- El as compared to the more massive SNe la. Of course this 
is the case if the changes are applied uniformly to the whole 
sample. In reality some events may have different values for 
the parameters listed in Eq when compared to each other. 

2 By plotting Amis (UVOIR) rather than the 56 Ni mass, we bypass 
the effect upon the luminosity (hence 56 Ni mass) associated with the 
uncertainty in the adopted distance to each event. 



The problem can, of course, be inverted to derive mean val- 
ues of q, k, and v e for a fixed ejected mass. With an ejected 
mass of 1.4 Mq we find mean values <v e > = 3762 km s _1 , 
<q> = 0.224, and <k> = 0.0080 cm 2 g -1 . If the two sublu- 
minous events (i.e. SN 1991bg and SN 1998de) are excluded, 
these parameters change to <v e > =3625 km s _1 , <q> = 
0.236, and <k> = 0.0084 cm 2 g" 1 . 

This e-folding velocity may be slightly on the high side 
compared to what is predicted from successful 1-D explosion 
models. However, it is not radically different from our adopted 
e-folding velocity. An explosion model with the majority of 
the 56 Ni mixed in the outer layers, as implied by a q = 0.224 is 
most unlikely for the vast majority of robust explosion models. 
But this could be the case for a Chandrasekhar mass progenitor 
that produces a subluminous SN la. As mentioned before, in 
the optically thin limit k ranges from 0.025 to 0.033 cm 2 g _1 . 
A factor of two (or even three) less of k — 0.025 cm 2 g _1 is 
unlikely. At most one could conceive of k to vary by ~ 50%. 

To calculate the amount of intermediate mass elements 
(IMEs) produced during nuclear burning, we simply subtract 
the ejected mass from the amount of 56 Ni produced. Note this 
this value also includes any stable Fe group elements and re- 
maining amount of unburned carbon and oxygen. These values 
are listed in Table 1. Excluding SN 1994D we find that the 
IMEs range from - 0.11 to - 0.67 Mq. 

In Fig.|5]we plot the ejected mass versus t , while holding 
the other parameters of Eq. @ constant. The solid line (case 
1) corresponds to all the fiducial values used to determine the 
ejected masses in Fig. [4] The dashed line (case 2) shows the ef- 
fect of keeping q and n fixed at the fiducial values while using 
v e = 3625 km s _1 . For the dash-dot-dot line (case 3) we used 
v e = 3625 km s -1 , K = 0.0084 cm 2 g" 1 and q = 0.5. Finally 
the dash-dot line (case 4) corresponds to v e = 3625 km s _1 , 
q = 1/3, and k = 0.0084 cm 2 g _1 . This figure illustrates the 
strong dependencies of the ejected masses. Masses much above 
the Chandrasekhar mass are achieved for only extreme cases. 
Both case 1 and case 2 provide ejected masses at or near the 
Chandrasekhar mass for events with large values of t and sub- 
stantially less for those events with values of t D R* 22-26 days. 

Another interesting feature displayed in Fig.|4]is that there 
appears to be little or no correlation between the ejected mass 
and the amount of 56 Ni. This is not entirely unexpected because 
as even with the presumption that all SNe la originate from a 
Chandrasekhar-size white dwarf, there still exists a range of 
ten or more in amount of 56 Ni produced. Nevertheless this is 
additional evidence which suggests that there is a significant 
variation in the burning of SNe la. 

We now turn our attention to the issue of the 7-ray es- 
cape fraction. In Fig. [6] we present the 7-ray escape fraction 
as a function of time (determined from Eq. 0) for five of the 
SNe la in our sample. As the ejecta of the supernova expands 
there is an increase in the 7-ray escape fraction. This can be 
attributed to the decrease in the column density, which is ac- 
companied with the expansion of the ejecta. Most of the curves 
in this figure are accompanied by a 'bump' between 20 and 
40 days past maximum light. Also included in Fig. [6] is the 
7-ray escape fraction calculated from W7. The agreement be- 
tween W7 and our calculated 7-ray escape fraction curves is 
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encouraging, considering that we are not adjusting any param- 
eters. 3 For the first three weeks after maximum light the 7-ray 
escape fraction from the UVOIR light curves is unreliable, as it 
is based on the assumption of r <C 1, which clearly is not the 
case fort < t D . 

From Fig-EJit is clear that the 7-ray escape fraction evolves 
faster in time for less luminous events. This is confirmed by 
Fig. [7] where we plot the 7-ray escape fraction at sixty days 
past bolometric maximum light versus Ami5 (UVOIR). At this 
epoch ~ 10% more 7 rays escape in the least luminous SNe la 
than in the brightest events. However, between 20 and 40 days 
past maximum light the differences are even more pronounced. 
This is expected, as at these epochs the UVOIR light curve has 
not yet reached its linear decline. In addition, the morphology 
of the secondary maximum can vary radically from SN to SN 
JSuntzefd2003llStritzingerl2005l) . This may then have a signif- 
icant effect on the evolution of the 7-ray escape fraction during 
these phases. 

5. Discussion 

With the stipulation that the UVOIR light curve reasonably 
traces the true bolometric flux, from within the period of soon 
after explosion to one hundred days past maximum light, we 
have been able to derive constraints on several important global 
parameters that relate directly to the progenitor systems of 
SNe la. The appeal of our approach is that, with relative ease 
and simple assumptions, we have used existing data to gain a 
deeper understanding on the origins of SNe la as well as pro- 
vide sorely needed constraints on current models. 

As previously mentioned, it is commonly believed that 
SNe la are the result of the thermonuclear disruption of a C-O 
white dwarf. The premise that thermonuclear combustion oc- 
curs at the Chandrasekhar limit was invoked to address the is- 
sue of homogeneity. However, today it is well established that 
SNe la are not true standa rd candles as once thought in the 
past (e.g. Leibundgut 2004). 4 Therefore we now must carefully 
scrutinize the data at hand in order to find plausible explana- 
tions that can account for the radical differences observed be- 
tween different SNe la. 

We first showed that the 56 Ni mass (hence luminosity) 
ranges from ~ 0.1 to ~ I.OMq. This confirms results pre- 



viously attained by 


several similar and independent meth- 


ods (Bowers et al. 1997; Cappellaro e 


:tal.ll997:Contardo et alJ 


2000; 




2002; Suntzeff 2003|. But this result is 



quite disheartening from earlier assertions that SNe la are 
standard candles. If all SNe la do indeed originate from a 
Chandrasekhar mass, an immediate question then is: What 
physical mechanism(s) can explain this range in luminosity? 



There has been considerable effort on the part of model- 
ers to address this question. Yet they have been met with little 
success to identify what parameter(s) can be tuned in order to 
account for a factor of ten in 56 Ni mass. Obvious candidates 

3 Note that we have assumed a rise time to bolometric maximum of 
19 days. 

4 Recently iKrisciunas et alJ 120041) have presented evidence that 
SNe la appear to be nearly 'standard candles' in the near-infrared. 



that may affect the production of 56 Ni are the initial param- 
eters prior to explosion, e.g. me tallicity, central density and 
ignition mechanism(s). Recently, Roo ke & Hillebrandtl ll2004 
have shown that the C-O ratio has essentially no effect on the 
amount of 56 Ni produced from burning to NSE. If prior to ex- 
plosion there is a significant amount of alpha elements within 
the white dwarf, one may reasonably expect the production of 
more stable isotopes, thus re ducing the amount of 56 Ni synthe- 
sized (Br achwitz et alJ2 000). Moreover it has been shown that 
changes in the central density of the white dwarf do influence 
the robustness of the explosion. Nevertheless, it is unrealistic 
that any one of these parameters, or even a combination of the 
three, can account for a factor of ten range in the 56 Ni mass. 
In reality these parameters affect the production of 56 Ni by no 
more than ~ 20%. 

More likely to influence the amou nt of 56 Ni synthesized is 
the e xplosion mechanism itself (see IStritzineer & Le ibundgut 
2005, and references within for a more detailed discussion). 
Currently the explosion mechanism and the subsequent evo- 
lution of the burning front is in open debate, and varies 
from a subsonic deflagration to a supersonic delayed deto- 
nation. Today the best Chandrasekhar mass models predict 
5b Ni masses that range between ~ 0.40 to 0.60 Mq. Due 
to computational limita tions the state -of-the-art 3-D deflagra- 
tion models jReinecke et alJl2 002albl) d o not produce copious 
amounts of 56 Ni JTravagho et alJ l2004). and have appreciable 
amounts of unburned carbon and oxygen left over in the in- 
ner ashes (Kozma et al. 2005). Th e delayed detonation models 
llKhokhlovlll99lt IWooslevlll990t IWooslev & WeaveJ Il994al: 
iHoflich & Khokhlovll99fj) ! on the other hand, can account for 
some of the more luminous events, but this class of models re- 
quires an additional free parameter. This parameter is essential 
to force the transition of the flame propagation from a deflagra- 
tion t o a detonation, and is physically not understood (how ever, 
see lGamezo et alJ2004l and lCiolombek & Niemeverl2005l) . The 
fact that there does not exist a single class of Chandrasekhar 
mass models that can account for the complete population of 
SNe la is quite dissatisfying and should be seriously addressed 
by theorists, if we are to insist that a Chandrasekhar-size white 
dwarf accounts for the progenitor system of all SNe la. 

Under the main assumption that at times greater than fifty 
days past maximum light the energy deposition in the ejecta 
of a SN la is solely due to the 56 Co— > 56 Fe decay chain, and 
thus the optical depth has a t~ 2 dependence, we can estimate 
(from the UVOIR light curve) the epoch when the photosphere 
transforms from being optically thick to optically thin. With 
th is know ledge we can then use the parameterized SN la model 
of lJeffervN 1999b to place constraints on the ejected mass. 

The results presented in Fig. [4] provide us with evidence 
that not all SNe la originate from a Chandrasekhar-size white 
dwarf or other very severe differences in the explosions like 
the distribution of 56 Ni or kinetic energies (expansion veloc- 
ities) exist. This would then immediately imply that some 
sort of sub-Chandrasekhar mass model is responsible for at 
least some SNe la. If true, this would be a radical change in 
thinking from the currently favored paradigm for the progen- 
itor systems of SNe la. However, the suggestion that a sub- 
Chandrasekhar mass model may be a viable candidate for the 
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progenitors of some SNe la is certainly not a new concept. 
Similar to Chandrasekhar mass models, previous attempts to 
simulate these systems have been plagued wit h their own prob- 
lem s. We refer t he reader to IHillebrandt & Niemeverl fcOOOh 
and Livi for detailed reviews concerning this class of 

progenitor system; we briefly summarize them here. 

Previo us attempts to model sub-Chandr asekha r ex- 
plosions dWooslev & Weaver! Il994bl Livne & Arnetil 1 19951 
iHoflich & Khokhlovll996r have met with some success in re- 
producing the observed light curves. However, these models 
typically predict a high-velocity layer of 56 Ni and helium above 
the IMEs, which is not observed in any spectra. It must be noted 
that relatively little effort has been made to conduct detailed 
3-D simulations of sMfr-Chan drasekhar mass models (but, see 
iGarcia-Senz etalJl 999 ; Benz 1997). With more detailed mod- 
eling, this progenitor channel may provide an attractive alter- 
native to the Chandrasekhar mass model. We also note that one 
appealing advantage offered by this model is the ability to ob- 
tain the progenitor statistics predicted by population synthesis 
calculations (s eelLiviol2000l and refere nces within). 

Previously, Cappellaro et al 1 11 19971) employed a technique 
that used observations of SNe la to determine both the 
56 Ni mass and the ejected mass. In their method they modeled 
the V^-band light curves of a small sample of SNe la using a 
simple Monte Carlo code. We find that our overall results are 
analogous to what they determined for both the range in the 
56 Ni mass and the ejected mass. Contrary to their work we em- 
ployed a different manner to determine these parameters and 
used UVOIR light curves rather than V-band light curves. By 
using the UVOIR light curve instead of the F-band light curve, 
we circumvented the crude assumption that the latter is a close 
surrogate to the former during post maximum times. Indeed, 
a comparison between our UVOIR light curves to the F-band 
light curves indicates that by fifty days past maximum light, 
the bolometric correction (mb r rn„) is ~ 0.2 mag or more. 
At later times this difference is amplified, as the near infrared 
passband s provide an increa sing contribution to the bolometric 
flux JSollerman et all2004l) . 

Although we find that our con clusions are in line with to 
those presented in Capp ellaro et all there are subtle differences 
between the four events that coincide in both studies. The num- 
bers we provide below for our results were obtained using 
Eq. @ and the fiducial values quoted previously. Also, note 
that there are slight differences (no larger tha n u = 0.20) in the 
distances used betwe en our work and that of fCappellar o et all 

For SN 1991bg. lCappellaro et"ai1 found a 56 Ni mass M Ni 
= 0.1 Mq and an ejected mass M C j = 0.7 Mq. This is quite 
comparable to our findings of M^i = 0.11 Mq and M C j = 
0.48 ± 0.14 Mq. Furthermore, our M Ni /M Gj ratio of 0.23 is 
larger compared to their 0.14. We found for SN 1992A, M Ni = 
0.40 Mq and M ej = 0.72 ± 0.27M Q , as compared to their 
M Ni = 0.4 Mq and M cj = 1.0 Mq. This then gives us a 
M Ni /M C j ratio of 0.56 compared to their 0.40. 

We find that our results with respect to the next two SNe la 
differ more than the first two stated events. For SN 1994D we 
calculated M Ni = 0.64 M and M cj = 0.65 ± 0.25 M , 
compared to their values of Mni = 0.8 Mq and M C j = 
1.4 Mq. Thus we obtain a larger difference in our MNi/M C j 



ratio of 0.98 compared to their 0.57. The fact that we have 
calculated a 56 Ni mass that is equal to the ejected mass is 
questionable. To determine th e 56 Ni mass of SN 19994D 
we used a new SBF distance jAiharetalJl200j^ rather than 
the S BF distance JTonrv et"ai1ll997|) used bv IContardo et alJ 
(2000) who determined a 56 Ni mass of 0.40 M ^. Using the 
distance modulus adopted by Capp ellar o et alJ we obtain a 
^ 6 Ni m ass of 0.67 MQ. Recently Feldmeie r. Jacobv & Phi llips 
( 2005) have calculated a planetary nebulae distance to the host 
galaxy of SN 1994D. In their study they have determined 
the distance modu lus /i = 30.66. This is comparable to the 
lTonrvetalJ ( ll997l) distances modulus /i = 30.68. Using the 
planetary nebulae distance the 56 Ni mass would be reduced 
to ~ 0.40Mq. None theless the 56 Ni mass determined by 
us and | C on|:ar(to et alJ is less than the 0.8 Mq calculated by 
Cappella ro et alJ with their method. The discrepancies between 
these values of the 56 Ni mass underscores the effect of the un- 
certainty in the distances. 

Finally, for SN 199 IT Cappell aroetail assumed M oj = 
MNi where MiMi = 1-1 Mq. We, on the other hand, found 
M Ni = 0.93 M and M cj = 1. 21 ± 0.36 Mq In summary 
we find the results presented by Cappellar o et alJ to be in fair 
agreement with our calculations, although some discrepancies 
do exist. 

We have presented an investigation of the bolometric be- 
havior of sixteen SNe la. In particular we have provided impor- 
tant constraints on the progenitor system(s) of these stellar ex- 
plosions. Our results suggest that some progenitor system(s) of 
SN la may emanate from the thermonuclear explosion of a sub 
Chandrasekhar-size white dwarf. This result may be difficult to 
reconcile with the current paradigm of the progenitor system 
of SNe la, i.e. a Chandrasekhar-size white dwarf. Moreover, 
our results suggest that the amount of 56 Ni produced during 
the explosion is most likely not dependent on the mass of the 
progenitor, but more likely on the manner in which nuclear 
burning is initiated and the subsequent dynamics of the flame 
propagation through the white dwarf. The range in synthesized 
56 Ni possibly indicates that there are two different explosion 
mechanisms. Further modeling of the explosion mechanism is 
required in order to investigate how different initial conditions 
can affect the observed range in luminosity. 

In Fig.|4]we see -in contrast to current thinking- that the 
mean ejected mass of many explosions is on the low side. 
A valid concern is that the parameters used to determine the 
ejected masses may not exactly represent those of a real SN la 
explosion. One parameter that could be in error, and does have 
a significant effect on our estimates of the ejected masses, is 
the adopted value of the e-folding velocity (see Fig[5j. If we 
assume slightly larger values of v e , the mean ejected mass for 
our sample would be in better agreement with 1.4 M©. How 
this parameter differs in 3-D simulations compared to 1-D sim- 
ulations is not yet clear. 

In addition, the simple assumption that any one of the pa- 
rameters in Eq. @ is unique for all events is probably incorrect. 
This may have a significant effect on the determined ejected 
mass for each event. However, this does not necessary imply 
that we would obtain larger ejected masses. It would be helpful 
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if the theorists in the future provided values of v e and q from 
their simulations. 

An acceptable argument concerning the results presented 
in this work is the validity of the model used to determine the 
ejected mass. There m ay be several a ssumptions built into the 
parameterized model o f l.TeffervHl999l) which may be too naive, 
and therefore the model may not adequately account for various 
complicated physical processes that occur within the progeni- 
tor of a SN la. However, there currently exists no other method 
to use observed photometry to place constraints on such a pa- 
rameter. 

In principle we would like to compare the UVOIR light 
curves to detailed NLTE modeled light curves. Unfortunately 
there has been little success in such an endeavor, owing to 
the complications in performing such time-dependent calcu- 
lations as well as t he limits of atomic line data, however see 
iKozma et all d2005l) . The next step will be to fit UVOIR light 
curves to a grid of model light curves produced from 3-D ra- 
diative transfer calculations, and then place further constraints 
on the progenitor systems of SNe la. 
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Table 1 . Well-observed SNe la 



SN 


Filters 


Ref. 


E(B-V)» al 


E(B-V) host 


/' 


vel 


Ref. 6 


M N i 


to 




Mime 














(kms- 1 ) 




(Mq) 


(days) 


CM©) 


CM©) 


SN1989B 


UBVRI 


1 


0.032 


0.340(0.04) 


30.22(0.12) 


797 


21 


0.64(0.18) 


32.23(0.12) 


1.06(0.32) 


0.42 


SN1991T 


UBVRI 


2 


0.022 


0.140(0.05) 


30.74(0.12) 


1012 


22 


0.93(0.30) 


34.44(0.23) 


1.21(0.36) 


0.28 


SN1991bg 


BVRI 


3, 4, 5 


0.040 


0.030(0.05) 


31.32(0.11) 


1322 


23 


0.11(0.03) 


21.62(0.11) 


0.48(0.14) 


0.37 


SN1992A 


UBVRI 


6 


0.017 


0.000(0.02) 


31.35(0.07) 


1341 


24 


0.40(0.03) 


26.58(0.10) 


0.72(0.22) 


0.32 


SN1994D 


UBVRI 


7, 8, 9 


0.022 


0.000(0.02) 


31.08(0.20) 


1184 


25 


0.64(0.13) 


25.31(0.11) 


0.65(0.20) 


0.01 


SN1994ae 


BVRI 


10 


0.031 


0.120(0.03) 


32.29(0.06) 


2067 


13 


0.84(0.13) 


32.33(0.13) 


1.07(0.32) 


0.23 


SN1995D 


BVRI 


10, 11 


0.058 


0.040(0.02) 


32.50(0.26) 


2272 


CMB 


0.66(0.23) 


35.15(0.12) 


1.26(0.38) 


0.60 


SN1995E 


BVRI 


10 


0.027 


0.740(0.03) 


33.42(0.18) 


3478 


CMB 


0.88(0.26) 


31.49(0.11) 


1.01(0.30) 


0.13 


SN1996X 


UBVRI 


10, 12 


0.069 


0.010(0.02) 


32.40(0.20) 


2174 


12 


0.73(0.21) 


28.70(0.10) 


0.84(0.25) 


0.11 


SN1998aq 


UBVRI 


13 


0.014 


0.002(0.05) 


31.66(0.09) 


1547 


13 


0.68(0.18) 


28.88(0.17) 


0.85(0.25) 


0.17 


SN1998de 


BVRI 


14 


0.060 


0.000(0.05) 


34.05(0.14) 


4653 


CMB 


0.09(0.03) 


27.80(0.10) 


0.68(0.20) 


0.59 


SN1999ac 


UBVRI 


15 


0.046 


0.120(0.05) 


33.06(0.21) 


2949 


CMB 


0.67(0.29) 


33.24(0.11) 


1.13(0.34) 


0.46 


SN1999dq 


UBVRI 


16 


0.024 


0.139(0.05) 


33.74(0.16) 


4029 


CMB 


0.80(0.29) 


34.91(0.14) 


1.24(0.37) 


0.44 


SN2000cx 


UBVRI 


16, 17, 18 


0.082 


0.000(0.05) 


31.90(0.20) 


1727 


23 


0.38(0.16) 


25.40(0.09) 


0.66(0.20) 


0.28 


SN2001el 


UBVRI 


19 


0.014 


0.206(0.05) 


30.83(0.54) 


1053 


CMB 


0.40(0.38) 


31.94(0.12) 


1.04(0.31) 


0.64 


SN2003du 


UBVRI 


20 


0.010 


0.000(0.05) 


32.23(0.30) 


2011 


CMB 


0.38(0.21) 


32.16(0.13) 


1.05(0.32) 


0.67 



a Taken from lSchlegel et alj Jl99^) dust maps. 

b For events without a direct distance estimate we selected Heliocentric velocities listed in NED and transformed these to the CMB reference frame. To 
account for peculiar velocities we assume throughout this work an error of 300 km s" 1 for all CMB distances. 
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n -i 

20 40 60 80 100 20 40 60 80 100 

days since L max 



Figure 1 Fit of Eq. (2) (solid line) to the UVOIR bolometric 
light curve (dotted curve) between 50 and 100 days past max- 
imum light. The dashed-dotted line is the energy deposition of 
7 rays and positrons from the 56 Ni to 56 Co to 56 Fe decay, as- 
suming complete trapping (i.e. r 3> 1). Dashed line is the case 
for complete escape of 7 rays (i.e. r <C 1). The vertical dotted 
line indicates the epoch (+50 days) when the fit begins. 

44 1 ' 1 1 1 1 1 ' ' 111 




40 h , , , , , , , , , , , , iii J 

20 40 60 80 100 

days since L 

Figure 2 Plot of Eq. (2) for a fixed 56 Ni mass of 0.38 M while 
varying t (solid lines). Here t ranges (from top to bottom) 
oo,45, 40, 35, 32.16, 25, 20, 15, and days. The light curve 
corresponds to SN 2003du. 
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Figure 3 Fiducial time, t plotted versus Ami 5 (UVOIR). Note 
the error bars associated with the values of t are smaller than 
the size of the points. 
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Figure 4 Ejected mass plotted vs. 56 Ni mass for 16 SNe la. 
Units are in solar mass. See text for comments concerning the 
error bars. 

Solid horizontal line indicates the Chandrasekhar mass. Slanted line 
has a slope of 1 . 
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Figure 5 Ejected mass plotted vs. t for fixed values of the 
parameters in Eq (0}. See Sect. 0]for a complete description 
of each curve. Solid vertical lines indicates the minimum and 
maximum values for t in our sample, and the solid horizontal 
lines indicates the Chandrasekhar mass. 
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Figure 6 7-ray escape fraction as a function of time since max- 
imum light for five SN la in our sample. These include (from 
top to bottom) SN 1991bg, SN 1994D, SN 2003du, SN 1991T, 
and SN 1999dq. The red dash line corresponds to W7. The ver- 
tical dot lines are the minimum and maximum values of t . For 
W7 we have assumed a rise time to bolometric maximum of 19 
days. 
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Figure 7 7-ray escape fraction at sixty days past maximum light 
vs. Am 15 (UVOIR). 



